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Abstract 

We discuss how momentum conservation is implemented in perturbative computa- 
tions based on expansions of anomalous dimensions appropriate at small x. We show 
that for any given choice of F% coefficient functions there always exists a factorization 
scheme where the gluon is defined in such a way that momentum is conserved at next 
to leading order. 
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An accurate determination of the x and Q 2 dependence of structure functions at small 
x requires solution of next to leading order (NLO) evolution equations appropriate to this 
kinematical region JIJ. The perturbative expansion of the anomalous dimensions or split- 
ting functions which govern perturbative evolution can then be reorganized in order to keep 
into account the presence of two large scales in the problem (Q 2 and s = (1 — x)Q 2 /x) @. 
The results of this procedure must be consistent with the constraints imposed by conser- 
vation laws. Energy-momentum conservation is particularly subtle in this context because 
energy is now one of the parameters which organize the perturbative expansion: therefore, 
it can only be imposed by an appropriate definition of the infinite set of contributions which 
are summed up in the anomalous dimensions. The problem of momentum conservation can 
thus be tackled only after these contributions have been defined in the most general way, 
in particular by correctly accounting for the freedom of choosing a factorization scheme. 
Here we will show that momentum conservation can always be implemented in NLO com- 
putations performed in expansion schemes appropriate to small x by a judicious choice of 
factorization scheme. After defining the appropriate small-x expansions, we will construct 
the most general transition functions which perform a scheme change. We will then show 
that given a matrix of anomalous dimensions and a factorization scheme, specified by the 
coefficient functions which relate F% to parton distributions (such as, for instance, MS or 
DIS) there is still enough freedom to perform a further scheme change within the given 
scheme (i.e. without changing the coefficient functions) such that momentum conservation 
is then obtained consistently at NLO. 

A determination of the evolution of parton distributions by solution of the renormal- 
ization group equations corresponds to summing leading (and subleading) logs of the form 
«s(log Q 2 ) 9 (log . In the usual loop expansion, appropriate to the Bjorken limit, the 
leading logs are those of Q 2 . Thus at leading order all terms with p = q are summed, 
at NLO those with q < p < 2q, and so forth; it then turns out that both at LO and 
NLO (and, in fact, at any order) < r < p. At small x, however, logs of - should also 
be considered leading, and the perturbative expansion reorganized accordingly. It is for 
example possible [[2] to define an expansion scheme appropriate to the Regge limit where 
the roles of In - and InQ 2 are interchanged, so that at LO all logs with p = r are summed, 
while 1 < q < p (the small x expansion) . A scheme where the two logs are treated on the 
same footing (double leading expansion) can also be defined, in which at LO any power 
of a s is accompanied by either of the two logs (or 1 < q < p, < r < p, 1 < p < q + r) , 
as well as a number of intermediate schemes. In the sequel we will consider specifically 
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the small x expansion, which is a theoretically interesting limiting case, and the double 
leading expansion, which is most interesting for phenomenology in the HERA region ||. 

The constraint imposed by momentum conservation on anomalous dimensions is a 
particular case of the general requirement that anomalous dimensions of conserved (or 
partially conserved) operators must vanish, which in turn is a consequence of the renor- 
malization group equations. This constraint takes the form 

y qg (l,a) + y gg (l,a) = 0, 7 W (1, a) + J gq (l, a) = 0, (1) 

where the anomalous dimensions are moments of the splittings functions, r )ij{N,a) = 
J Q x n Pij(x;t), and depend on t = ln(Q 2 /A 2 ) through the running of the coupling a. 
Momentum conservation [eq. (jl])] is imposed order by order in the usual loop expansion in 
powers of a of the anomalous dimensions, by suitable choice of normalization of the quark 
and gluon distributions. If, however, the order a k contribution to 7^ is further expanded 
in powers of iV the single terms of this expansion will not, of course, satisfy eq. (]!]). Now, 
the various expansion schemes alluded above are obtained precisely by performing such 
expansions, and then including at each order a suitable subset of terms, which then will 
not conserve momentum automatically. 

When evolution in the small x region is approached by choosing an appropriate ex- 
pansion scheme, momentum conservation should be imposed order by order, just as it is 
in the loop expansion.^ That this is a priori non-trivial is clear from the observation that 
the leading order anomalous dimensions in the small x expansion described above violate 
eq. ([[]). Just like in the usual expansion, however, eq. (|l|) only holds if parton distributions 
are defined appropriately, i.e. for suitable choices of the factorization scheme. Before we 
discuss the implementation of momentum conservation we must therefore discuss changes 
of factorization scheme within various small x expansions. 

Changing the factorization scheme amounts [| to a redefinition of the singlet par- 
ton densities f(N,t) = (^), where q = YLiiQi + Letting /—>/' = Uf, the naive 
partonic interpretation of the parton densities will be maintained if we always assume 
that U(N,a) = 1 + 0(a). The renormalization group equation 4rf = 7/ then remains 
unchanged provided 7 — > 7', where 

i = U 1 U~ 1 + (f t U)U- 1 . (2) 

1 Alternatively, momentum conservation could always be imposed at each order in a by in- 
cluding an Ansatz for the (yet unknown) terms which are formally sub-subleading in the small x 
expansions: several proposals of this kind are discussed in ref. S- 
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Because U only depends on t through a 

^U(N,a)=P(a)^U(N,a), (3) 

where (3(a) is the beta function. Since (3(a) = —(3oa 2 + 0(a 3 ) it follows that the second 
term of (0) is subleading compared to the first. 

We must now choose a specific expansion scheme. We will eventually prove our result 
in the physically relevant double leading expansion; however, we consider first the small 
x expansion, since this will allow us to present the general structure of our results in a 
somewhat simpler setting. In the small- x scheme the singlet anomalous dimension are 
given by 

7 = 7 S + 7 SS + ---; 

oo 

7 s (iV,a)^ 7>/^) n 

n=no V / 

oo 

Y s (N,a) = a J] 7n S +i(«/A0 n , 

n=no — 1 

where 7 s (N, a) sums the leading singularities, 7 ss (iV, a) the subleading singularities, and 
so on. The one loop contributions to 7 s and 7 SS [i.e. the contributions with n = 1 in 
eq. @] will in general violate the condition ([!]); they are scheme independent and thus 
lead inevitably to momentum nonconservation in the usual small x expansion if no = 1. 
We will hence consider here a truncated small x expansion with no = 2, i.e. with these one 
loop contributions suppressed. The results then found will coincide with those of the usual 
small x scheme in the x — > limit; we will use them to prove that momentum conservation 
holds in the double leading expansion. 

At NLO in a typical factorization scheme such as the QoDIS scheme [g| (or the MS 
and DIS schemes @) 7 is then of the form 

1= ( n 9 7,W W n? 7fV ... (5) 
\n 9 + i igj \n s g igj \n s g s + i 9 i 8 9 s J 

where r = Cf/Ca is the colour-charge factor. It is important for what follows to notice 
that both of the quark anomalous dimensions vanish at leading order; also, both the gluon 
anomalous dimensions at leading order, and the quark ones at subleading order obey a 
colour-charge relation.! At subleadina; order in the gluon sector the anomalous dimensions 



2 Notice that this relation is however not satisfied by the one-loop term "y( s , which has been 
suppressed. 
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7|« and 7!* are as yet unknown beyond two loops so we introduce 7 g = 7^ — r^ s g s g to allow 
for possible violations of the colour-charge relation at NLO in the gluon channel. 

Consider now parton distributions normalized so as to satisfy the momentum sum rule 
at some scale t = 0: q(l, 0) + g(l, 0) = 1. Momentum is then conserved in the evolution of 
these distributions to the scale t if the two conditions ([]]) are satisfied. With anomalous 
dimension of the form ([5]), the momentum sum rule will be violated at LO (and thus in a 
scheme invariant way) since 7^(1, ot) 7^ 0. However, the quark distribution only evolves at 
NLO, and the gluon can only be directly observed at NLO (through measurement of Fl, 
say); therefore, this LO violation has no physical effects. Momentum conservation starts 
thus being physically relevant at NLO, where it imposes a condition relating the LO and 
NLO components of the anomalous dimension, as well as the colour-charge relation in the 
gluonic sector: 

(7 9 SS )n + (lg)n + (7f )n = 0, (%) n = 0, (6) 

where 7^ = J2(lg)n(&/N) n , etc. The conditions (|6|) will not in general be satisfied in a 
given generic scheme; however in what follows we will show that it is always possible to 
find factorization schemes such that both conditions are satisfied (for n > 1), and thus in 
which momentum is conserved. 

Consider first a scheme change U which is LO in the small- a; expansion, i.e. U = 
1 + U n (a/N) n . The most general form of U which retains the identification of F<i with 
the quark density at LO in this expansion is§ 

U=(l U ), (7) 
\u u J 

where u = 1 + J2T U n((x/N) n while u = J2T ^n(a/iV) n and similarly for u. Substitution 
in (0) gives to LO 

/ / f u{ru - u) u(l — ru)\ n/ , , . 

7 =^ g {u-uu) ) ' ) ' )+0(a). (8) 

y yu[ru — u) u(l — ru) J 

If we now insist that the leading order anomalous dimension is to remain unchanged, we 
must choose 

u = 0, u = r(u-l). (9) 



3 We neglect LO transformations proportional to the unit matrix, since they modify the LO 
relation between F2 and q. 
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A LO scheme change thus amounts essentially to a redefinition of the gluon normalization 
by the function it. The NLO anomalous dimensions are then correspondingly modified 
according to 



(10) 



so that in particular r y qg = ^ q ^> 'Jq/u. Notice that the colour-charge relation is automati- 
cally preserved not only in the LO gluon sector, but also in the NLO quark sector. 

The quark distribution is left unaffected by the LO transformations considered so far. 
However, we may still perform a NLO redefinition /'—>/"= (1 + V)f, where 

oo 

V = aJ2Vn+i(a/N) n . (11) 
o 

This induces a corresponding change in the anomalous dimension 

i^i' = l'+[V n '] + 0{a 2 ). (12) 

Writing 

V=(* V ), (13) 
and keeping only LO and NLO terms, we find 

7" = 7' + 7i(, ^ ~, V )+0(a 2 ). (14) 
y \[rw — w — rv) —rv J 

Combining this NLO transformation with the LO transformation U, we thus have alto- 
gether at NLO 

Iqq = 111 + r ( U ^ ~ !)7q S + ™7p, 

iq 9 = 7« + - ihr + 

Igq = 7 3g + (u- + r%) - r 2 (w _1 - l)7^ s + (rw - w - 7-6)7* + rf t lnu, 
I99 = ^99 - riu- 1 - l)Y q s ~ rvY g + ft lnu - 

The colour-charge relation in the quark sector is thus again preserved automatically. 
Starting for definiteness in a parton scheme (such as, say the QoDIS scheme), where 
F 2 (x,t) = (e 2 )xq(x,t) (and (e 2 ) is a numerical factor) the functions v and v then give the 



F<i coefficient functions, since q = (1 — v )q" — vg" . To enforce the coiour-charge relation 
in the coefficient functions we should take 

v = rv + v; v = —Via (16) 

(as is the case in MS scheme, for example Q). Similarly, the NLO terms violating the 
colour-charge relation at NLO in the gluon sector can be removed by choosing rw — w 
appropriately 

The most general NLO scheme change is thus parameterized by two parameters u 
and v (or more properly a two-fold infinity of parameters, u n and v n ), plus a parameter 
v which is fixed requiring the coefficient functions to satisfy the colour charge relation, 
and a parameter rw — w which can be used to impose the colour-charge relation in NLO 
anomalous dimensions, while the orthogonal combination w + rw has no effect at all at 
NLO. Given then parton distributions in a particular parton scheme (such as QoDIS) 
the parameter u redefines the normalization of the gluon, without affecting F2 directly, 
thus takes to different parton schemes (such as DIS or SDIS) while v moves to a different 
(generally non partonic) scheme (such as MS). 

We can now extend these results to the physically relevant double leading expansion, 
which treats the two large scales on an equal footing (other expansion schemes can be 
handled in a similar way). In this expansion, the matrix of singlet anomalous dimensions 
consists of a large x and a small x contribution, j L and 7 5 , respectively: 

7 = 7 L + 7 S ; 7 L = 7'+7 2 , 7 S = 7 S + 7 SS , (17) 

where 7 1 and 7 2 are the usual one and two loop anomalous dimensions, while the leading 
and subleading singularities have the form (|4|), but now with no = 3, so that both the 
one and the two loop terms are suppressed to avoid double counting. Scheme changes at 
LO are still effected by the matrix U eq. ([7]) (no nontrivial LO scheme change of 7 L is 
possible). These transform the singular anomalous dimensions according to eq. @, (|T0| ) 
and are thus subject to the constraint eq. (^). In addition, further singular contributions 
are produced by the action of U on 7 L , i.e. U"f L U~ 1 . These may be combined with the 
LO transformation eq. ( |I0D of 7 s to give 

l' qg = l qa + - l)7w, 

Igq = Igq + ( u ~ 1) [(lm - r lfgg + r lqq) + ^ ~ ^)lqg] + r jt lliU ' 

%g = ^gg ~ r^' 1 - l) 7w + | lnu, 
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independently of the split of 7 into r ) L and 7 s . Notice that all terms on the right hand side 
beyond the first are NLO (because in particular wy 1 = 0(a(a/N) n )), and that neither r y 2 
nor 7 s make any contribution to them. 

Again the LO transformation redefines the normalization of the gluon distribution 
without affecting the quark, which only transforms under a NLO transformations. This 
now contains two components: a small x NLO transformation of the form (|TTD, fll3| ) 
considered previously, plus a standard scheme changing 0(a) function V L , which vanishes 
as N -> 0: 

00 

V = V L + V S ; V s = aY,V: +1 {a/N) n . (19) 

n=0 

Since V L has no effect on 7 s at NLO, producing only terms of 0(a 2 (a/N) n ), we will not 
consider it any further here: its only effect is to change the two loop anomalous dimension 
7 2 in the usual way. The effect of the remaining singular transformation V s is then very 
similar to that already discussed in the small x expansion: 7 s transforms according to the 
small x eq. (|I4[), while the 0(a) contribution to V s (which could in fact equivalently be 
viewed as a contribution to V L ) also produces a NLO change of 7 L . The general scheme 
transformation at NLO in the double leading expansion is thus 

Y' = y + Y g ( ( ^ ~x V )+a[V 1 s , 1 1 }, (20) 
9 V (rw — w — rv) —rv J 

where here 7* is given by (|J) with no = 1 (i.e. including the one-loop contribution) v was 
defined in eq. flT6|), and 7' is the anomalous dimension transformed at LO according to 
eq. ®. 

We thus find again that, starting with a parton scheme, the parameter v takes us to 
nonpartonic schemes, while the further parameters rw — w and v may be fixed by requiring 
the color-charge relations in the anomalous dimensions 7^ and in the corresponding small 
x coefficient functions, respectively. In comparison to standard scheme change at large 
x Q there seems thus to be an additional freedom, parametrized by u, of redefining 
the normalization of the gluon. This redefinition by a LO function of a/N modifies the 
singular part of the NLO anomalous dimension (and appears thus to be peculiar of small x 
expansion schemes). The effect of this transformation can in practice be rather important. 
Two (in some sense extremal) choices have been considered in the literature. Taking 
u(N, a) = R(N, a)^ 1 , the singular normalization factor for the gluon in MS ||, takes one 
from the QoDIS scheme to the more singular DIS scheme 0. Because of the singular 
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behavior of R(N, a) the anomalous dimensions 7^ s in this scheme are substantially larger. 
Conversely, one may take u = ~fqg/~f\ g - this removes the singular terms in the quark sector 
altogether, factoring them in the initial gluon distribution (the 'SDIS' scheme^ [[7]). 

The effect on the leading order relation between Fl and the gluon distribution is 
equally significant. In the double leading expansion and QoDIS factorization scheme || 
Fl = (C L + liL)g + 0(a 2 ), where C\ is the usual two loop longitudinal coefficient function, 
and hi, = ceJ2T hn+i(oi / N) n ; after the (LO) change of scheme 

F L = (Cl + fti)(«-V + riu- 1 - l)q') + 0(a 2 ). (21) 

The result is especially clear in the small x expansion, where the quark distribution is 
subleading compared to the gluon 0: the dominant contribution to Fl is then equal 
to F L = h L Rg' + 0(a 2 ) in DIS @, while in SDIS F L = (h L /h 2 )c q ag' + 0(a 2 ) (since 
7* = c q a + 0(a 2 )). Because (h L /h 2 ) = (1 - 7 *)/(l + § 7 *(1 - 7 p) = 1 + 0(a/N), at 
small x Fl will then differ from the gluon distribution by a large factor. The leading order 
identification of Fl with the gluon may however be restored by choosing a 'GDIS' scheme 
in which u = 1 + Hl/C^ so that Fl = C\a.g' , so that in the small x expansion Fl = 
c q ag' + 0(a 2 ) since C\ = c q a + 0(a 2 )). In the GDIS scheme 7 g = (h2/hL)c q a + 0(a 2 ), so 
the quark anomalous dimensions are less singular than in QoDIS or DIS, but more singular 
than in SDIS. 

This extra freedom in the choice of factorization scheme at small x can however be 
constrained by requiring momentum conservation. Combining the two transformations 
( |18D and ([2(]) and imposing the constraint of momentum conservation ([!]) we get at N = 1 
the two conditions 

-Poa 2 £u + (Y g + 7f + n s q s - (r - l)v%)u = (r - l)<yj', (22) 
7 S + (u - 1)% + (rw -w - rv)jg = 0, (23) 

where we used the crucial fact that momentum conservation at one and two loops, i.e. 
7 1 (l,a) = 7 2 (l,o;) = naturally eliminates all non-singular terms. Notice that Eq. (|23| ) 
is equivalent to the color-charge relation (in the subleading gluon sector) , which then does 
not have to be imposed separately, but rather follows automatically. 



4 One could also consider a SDIS' scheme @, where u = 7qg/(7g g + 7q g ), so that j q reduces 
to its two-loop expression 
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Now, for each choice of u eq. (22) has a unique solution for v[u] [i.e. for the infinite 
set of coefficients v n (u m )], while similarly ( ]2"B| ) has a unique solution for (rw — w)[u\. Con- 
versely, since for a given v eq. (p2|) is a first order differential equation for the function 
u(l, a) with boundary condition u(l,0) = 1 (compatible with ( p2|) because by construc- 
tion 7 s and 7 SS contain no one loop terms), it always has a unique solution u[v], so the 
functional v[u] is invertiblei! It follows that eq. (p2|) defines a monotonic curve u[v] in the 
two dimensional space of schemes (u, v) along which momentum is conserved: imposing 
momentum conservation at NLO the gluon normalization is fixed uniquely by the choice 
of F2 coefficient functions.il 

In order to actually determine the curve u[v], however, one needs full knowledge of the 
(unknown) subleading gluon anomalous dimension, i.e. of 7* s and 7 3 in a given scheme, 



say QqDIS. Indeed, the same argument which shows that for a given 7 S (l,a) eq. (|2 



determines v[u] or u[v] also implies that for a given v eq. (|22| ) determines it[7j s ] (or 7g S [u]): 
for every v and 7^ s there exists a u which conserves momentum. Conversely, for every pair 
of u and v momentum conservation determines a unique two dimensional surface 7| s [w, v], 
whose intersection with the plane 7* s = (7g S )Q DIS §i ves back the curve u[v]. Moving on 
the plane parametrized by u and v for fixed 7^ yields anomalous dimensions which are 
equivalent up to a NLO change of scheme (and only conserve momentum along the u[v] 
curve); thus, it explores the uncertainty related to the ignorance of NNLO corrections. 
Different choices of 7^ s then explore the uncertainty related to the ignorance of the NLO 
gluon anomalous dimension. 



5 While always true in QCD with iV c colours, since then r = \{1 — 1/N%), in supersymmetric 
Yang-Mills r = 1, so ( p2|) fixes u independently of v. 

6 One may, of course, leave the line v[u] while preserving momentum conservation at NNLO: 
for example one could perform a scheme change of the conventional momentum conserving form 
V = (-™ 7) (where v is NLO). Under this change of scheme 7 — ► 7', where, if 7 satisfies the 
colour-charge relation in both quark and gluon sectors 

7' = 7-£(7 9 +7,)(_ r r 2 _M+0(a 2 ). (24) 

If 7 satisfies momentum conservation, (7^ + 7 9 )jv=i = and the second term vanishes; however 
vy q , which must be kept in order to get momentum conservation, is actually NNLO. One would 
then be forced to include some NNLO contributions in a NLO computations in order to preserve 
the momentum sum rule. 



9 



Momentum conservation constrains the uncertainty in that it fixes one of the param- 
eters in terms of the other two; the overall uncertainty however is still necessarily NLO. 
Now, given Y g s and r y' q ss (from eq. ( [18D for a certain choice of u and v ) Y g ss is determined 
algebraically (by eq.(|l|) for n > 2) without the need for an explicit computation ||; it is 
thus convenient to fix it thus, and vary u for given v in order to estimate the corresponding 
NLO uncertainty [|l|. In practice, this uncertainty should only have minor effects on F2, 
because since 7* vanishes, q evolves only at NLO, while although g evolves at LO it only 
affects structure functions at NLO; the effect of 7| s is then subleading when compared to 
7!, and thus effectively NNLO. 

In conclusion, we have shown that momentum conservation can be enforced in QCD 
evolution equations at NLO even when anomalous dimensions are computed within an ex- 
pansion scheme which sums up all leading and subleading logs of both ^ and Q 2 , such as 
the double leading expansion 0] . Within these expansion schemes there is a wider freedom 
of choice of factorization scheme than in the usual (large x) loop expansion of anomalous 
dimensions J|] : besides the usual freedom to perform a NLO scheme change which modifies 
the F2 coefficient functions, there is now also the possibility of performing a LO scheme 
change which does not affect the F2 coefficient functions or the LO anomalous dimensions, 
but changes the definition of the gluon distribution (or, equivalently, the Fj_, coefficient 
function). The latter freedom is however fixed by the requirement of momentum conser- 
vation. This is analogous to what happens in large x computations, where the momentum 
sum rule fixes the normalization of the gluon; however at small x the normalization is actu- 
ally given by a general LO function of a/N, rather than being just a number. Because the 
NLO 7" and 7 s9 gluon anomalous dimensions are still unknown, this freedom cannot in 
practice be pinned down, so that it is still necessary to vary the gluon normalization within 
a given NLO factorization scheme, thereby introducing a NLO uncertainty in the solution 
of the evolution equations 0] . This uncertainty has fortunately only sub-subleading effects 
on the structure function F%, however, due to the vanishing of the LO singularities in the 
quark anomalous dimensions. 

If 7| s were to be determined explicitly (for instance by computing subleading correc- 
tions to the BFKL kernel in a particular factorization scheme) it would then be possible use 
momentum conservation to fix the definition of the gluon distribution, thereby reducing the 
uncertainty in the solution of evolution equations to a purely NNLO one. This could have 
significant phenomenological consequences in that it might substantially reduce the large 
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scheme dependence which is at present an intrinsic feature of perturbative computations 
at small x and relatively low Q 2 ||. 

Acknowledgement: We would like to thank Keith Ellis for emphasising to us the 
importance of momentum conservation at small x. 
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